SENSOR FOR DETERMINING THE ANGULAR POSITION 
OP A RADIATING POINT SOURCE IN TWO DIMENSIONS 
AND METHOD OF OPERATION 

BACKGROUND OF THE INVENTION 

The present invention relates to a sensor for determining the 
angular position of a radiating point source in two dimensions and its 
method of operation. The present invention has particular pertinence in 
the field of surgery, with particular reference to close range surgical 
instrument tracking. Additionally, the present invention is applicable in 
the fields of virtual reality and pilot head tracking systems. 

In the prior art, systems have been devised that employ one 
dimensional masked or coded apertures placed above multi -element one 
dimensional detectors. Such structure is disclosed in U.S. Patent Nos . 
4,810,870 and 5,408,323. In such systems, correlation techniques are 
employed to determine the position of a projected pseudo-random pattern on 
a detector surface by comparing signals from fixed detector elements 
arranged in a unique pattern or from comparison with fixed reference 
signals previously stored in a computer. 

Additionally, it is known in the prior art to use a two dimensional 
aperture and a one dimensional detector in the form of a charge coupled 
device (CCD) to determine angular position in two dimensions. Such a 
system is disclosed in U.S. Patent No. 4,092,072. A drawback of such a 
system is that it does not provide immunity from errors resulting from 
reflections, and does not yield sufficiently high enough resolution as is 
the case with techniques employing pseudo-random patterns and correlation 
techniques . 

The prior art systems described above, while capable of determining 
angle of incidence with high resolution, yield only one -dimensional 
information. The present invention differs from the teachings of these 
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patents as providing high resolution two-dimensional information. It is 
with this difference in mind that the present invention was developed. 

SUMMARY OF THE INVENTION 
The present invention relates to a sensor for determining the 
angular position of a radiating point source in two dimensions and its 
method of operation. The present invention includes the following 
interrelated objects, aspects and features: 

(1) In a first aspect of the present invention, the present 
invention includes a method for determining the angular position in two 
dimensions of a point source of radiation with respect to a detector by 
examining shifts of periodic components, in two dimensions, of a projected 
image of a variable transmissivity mask. In the preferred embodiment of 
the present invention, the mask is encoded in two skewed directions with 
waveforms consisting of several frequencies in prescribed patterns. 

(2) The frequency spectra of the received detector patterns are 
computed. In order to facilitate such computations, the constituent 
frequencies are substantially separated so as to be distinguished in the 
Fast Fourier Transform (FFT) . 

(3) Each of the frequency patterns that are coded on the variable 
transmissivity mask consists of a series of low frequencies followed by a 
series of variable frequencies, and a series of high frequencies. 
Applicants have found that accurate measurements may be obtained through 
creation of FFT plots for various pattern shifts. The middle variable 
frequencies exhibit frequency changes responsive to various image 
positions. The low and high frequencies that straddle the variable 
frequencies have been found to be responsive in phase to variations in 
image position. 



(4) Applicants have found that the frequency variations in the 
variable frequencies may be used to indicate coarse position in each 
coordinate axis while the phases of the fixed low and high frequencies may 
be used to indicate medium and fine position in each coordinate axis. 

As such, it is a first object of the present invention to provide a 
sensor for determining the angular position of a radiation point source in 
two dimensions and the method thereof. 

It is a further object of the present invention to provide such a 
sensor in which a variable transmissivity mask is created by encoding the 
mask with sequences of frequencies in two skewed or orthogonal directions. 

It is a still further object of the present invention to provide 
such frequency patterns including use of low frequencies followed by 
variable frequencies followed by high frequencies. 

It is a still further object of the present invention to provide 
such a sensor in which coarse measurements of position may be carried out 
through measurements of frequency changes in the variable frequencies and 
medium and fine position measurements may be carried out through 
variations in the phases of the low and high frequencies. 

It is a yet further object of the present invention to provide such 
a sensor, in a second embodiment thereof, wherein the mask pattern is 
formed by a first function including low, variable and high frequency 
components, a second function with fixed low and high frequency 
components, and a third function with variable frequency components. 

These and other objects, aspects and features of the present 
invention will be better understood from the following detailed 
description of the preferred embodiments when read in conjunction with the 
appended drawing figures. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a schematic representation of a preferred embodiment 
of the present invention. 

Figure 2 shows a front view of the embodiment of Figure 1 showing 
the X and Y axes . 

Figure 3 shows a view rotated 90 degrees from the view of Figure 2 
to show the X and Z axes. 

Figure 4 shows a top view of the embodiment of Figure 1 . 

Figures 5, 6 and 7 show graphs of computed Fast Fourier Transforms 
(FFTs) for first, second and third pattern shifts, respectively. 

Figure 8 shows a graph of a computed FFT showing both magnitude and 

phase. 

Figure 9 shows a schematic representation of the method steps of the 
present invention . 

Figure 10 shows a schematic representation of an alternative 
embodiment in which the mask pattern is placed in a configuration 
different from that of the embodiment best illustrated in Figure 1. 

SPECIFIC DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Reference is first made to Figure 1 which shows a schematic 
representation of a first embodiment of the present invention. In Figure 
1, a point source of radiation 1 is directed toward a mask 3. The 
radiation travels through the mask 3 and, as modified by the mask 3, is 
received by the detector 5 . 

With further reference to Figure 1, the mask 3 has a pattern thereon 
which is created by a first set of frequencies along the double-headed 
arrow 7 and a second set of frequencies along the double-headed arrow 9. 

Figures 2, 3 and 4 show front, side and top views, respectively, of 
the system illustrated in Figure 1. Figures 5, 6 and 7 provide a better 



understanding of the frequency sequences shown at the reference numerals 
7 and 9. Each of the frequency sequences consist of at least two low 
fixed frequencies followed by at least two variable frequencies followed 
by at least two fixed high frequencies. 

The detector 5 may, if desired, consist of a linear CCD detector 
with the mask 3 situated at a prescribed fixed distance therefrom and, as 
seen in Figures 2 and 3, the mask 3 is preferably in a plane parallel to 
the plane of the CCD detector array 5. As should be understood from 
Figures 1-4, the point source emitter 1 illuminates the mask 3 casting an 
image onto the plane of the detector 5. The detected image is recorded in 
a computer (not shown) and the Fast Fourier Transform (FFT) is computed. 

Applicants have found that variations in the frequencies of the 
variable frequencies shown in the middle of the graphs of Figures 5, 6 and 
7 indicate coarse position of the point source 1, while variations in the 
phases of the fixed low and high frequencies indicate medium and fine 
position. 

Since the point source is not located at a distance of infinity, the 
detected image pattern on the detector 5 is never equal to the dimensions 
of the image of the mask 3, but is always magnified with the magnification 
being as much as 1.05 to 1.10. The magnification factor 1^ is computed by 
taking the ratio of measured detector frequency to actual mask frequency 
preprogrammed into the computer. 

* The manner of measuring position in two dimensions in accordance 
with the teachings of the present invention will now be described in 
greater detail below. 

In the preferred embodiment of the present invention, the mask is 
programmed with two series of waveforms that are related to one another as 
being located in two axes skewed with respect to each other and with 
respect to the axis of the detector 5. Each waveform preferably consists 



of a plurality of fixed low frequency waves followed by a plurality of 
variable frequency waves followed by a plurality of fixed high frequency 
waves. In the example of Figures 1-4 in which the detector 5 lies on the 
y axis (see, in particular, Figure 4), one waveform is located along a 

5 line defined where y = z and the other waveform is located on a line 

defined as y = -z. These lines are orthogonal with respect to one 
another. The line defined by y =z may be described by the function g(y+z) 
and the line defined by y = -z may be described by the function h(y-z) . 
The transmissivity equation is: 

10 Mask(y,z) = g(y+z) + h(y-z) (D 

Replacing g(y+z) and h(y-z) by their constituent functions, the mask 
transmissivity is described by: 

2) a 05 + 1- ~ ^ 



Mask(y, 



(2) 



15 Where F lowl/ F low2 , F variablel , F variable2 , F highl , F high2 , are low, variable and high 

frequency constituents, respectively, of the mask, and modi, and mod2 are 
modulation factors. The detector lies on the y axis and has 1,024 
elements or pixels, making N=l,024 in the equation. 

Figures 5, 6 and 7 show FFTs for various pattern shifts with the 

20 middle two frequencies comprising variable transmissivity frequencies and 

being responsive in frequency to various image positions. The left two 
frequencies are fixed low frequencies and the right two frequencies are 
fixed high frequencies. These fixed frequencies are responsive in phase 
to various image positions. 

25 From FFT plots, the spectral peaks are determined. The position in 

frequency of the variable frequencies indicates coarse position in each 
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coordinate axis, while the phases of the fixed frequencies indicate medium 
and fine position in each coordinate axis. Figure 8 shows the FFT 
magnitude and phase of a fixed frequency component. 

Because the point source is not located at infinity, the detected 

5 image pattern is never equal to the mask image, but is always magnified, 

with magnifications as much as 1.05 to 1.10. The magnification (KJ is 
computed by taking the ratio of the measured detector frequency to the 
actual mask frequency. 

= Fhighl M8k /Fhighlae t = Fhigh2 mask /Fhigh2 det (3) 

10 The following image shift calculation uses a mask with variable frequency 

ranging from Fmin to Fmax over a range of 3N pixels. The coarse pixel 
shift indicated by the variable frequency component is: 

Y m _van = [Fvariablelaet*!^ - Fmidl^J / [Fmaxl niaBk - Fminl^] 

* [3*N] pixels (4) 

15 Y ravar2 = [Fvariable2 det *K m - Fmid2 mask ] / [Fmax2 mask - Fmin2 raask ] 

* [3*N] pixels (5) 

To compute the shift of the low and high frequency components, 
corresponding to medium and fine shifts, the phase at the array midpoint 
must be computed. The FFT routine returns the phase of frequencies at 0°, 
20 so to compute the midpoint phase, (i = 511.5) the following expression is 



used: 

Arg(Flowl det ) | 1-5lliS = Arg (Flowl det ) | Lo0 + Mod 2n (511 . 5*Flowl det *2n) (6) 
Arg(Fhighl det ) | io511 . 5 = Arg (Fhighl det ) | i-0 + Mod 2n (511 . 5*Fhighl det *2n) (7) 
25 Arg(Flow2 det ) | i=511 . 5 = Arg (Flow2 det ) | io0 + Mpd 2n (511 . 5*Flow2 det *2n) (8) 

Arg(Fhigh2 det ) \ Lm511 . 5 = Arg (Fhigh2 det ) | 1-0 + Mod 2n (511 . 5*Fhigh2 det *2n) (9) 

The shifts relative to mask functions g(y+z) and h(y-z) are indicated by 
the fixed frequency components by: 

Yin lowl = t^lowl 

(Arg(Flowl det ) |i. 511 . 5 /2n)] /Flowl^ (10) 

30 Y m hlghl = [P highl + (Arg(Fhighl det ) | i=sll . 5 /2n) ] /FhighW (11) 
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Y m _iov 2 = [Piow 2 + (Arg(Flow2 det ) | ia511 . 5 /2n)] /Flow2 mas)c (12) 
Y ra _ high2 = [Phigh + (Arg(Fhigh2 det ) | i=511 . 5 /2n) ] /Fhigh2 ma8k (13) 

Z m _highl - Y m highl (14) 
Z m _high2 = ~Y m ^ high2 (15) 

Where P lowl and P low2 are determined from the variable frequency shift 
estimations Y m varl and Y mvar2 , Phighi and P high2 are determined from the low 
frequency shift estimations Y mlowl and Y rolow2 . 

Piowi = Nearest_Integer [ (Y^ar^Flowl^sk) - 

(Arg(Flowl dec ) | i. 511 . 5 /2n) ] (16) 
P high2 = Nearest_lnteger [ (Y ni _ low2 *Fhigh2 mask ) - 

(Arg(Fhigh2 det ) | li5U . s /2n) ] (17) 

The overall coordinates of incidence to the mask are computed by: 

Y m = Y m highl + Y m _hi g h2 ( 1 8 ) 

Z ra = Z m _ highl + Z m _ high 2 

The angles of incidence to the mask are computed using: 

9^ = Tan" 1 (Y m /X m ) (20) 

e zm = Tan-Mz./x.) (21) 
Figure 9 shows a block diagram of the calculation process. 

As seen in Figure 9, the first step in performing the method of the 
present invention is to record the detector array image and, therefrom, to 
compute the image Fast Fourier Transform (FFT) . Next, the highest 
frequency F d is identified and then the magnification constant is 
computed. 

Coarse measurements are then made by determining the peaks of the 
variable frequency waveforms F variable and, thereafter, the phases of the low 
and high frequency waveforms are determined. From those determinations, 
total pattern shifts as a sum of the frequency component phase shifts are 
determined and, therefrom, pattern shifts of the magnified image are 
calculated. 
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From these measurements and calculations, the actual unmagnified 
mask pattern shift is computed and, with this computation having been 
completed, the angles of incidence for the point source of light are 
accurately computed. 

A second embodiment of the present invention is schematically 
illustrated in Figure 10. In Figure 10, the system 10 is seen to include 
a point source of radiation 11, a mask 3, and a detector 5. The mask 3 
has a pattern thereon including fixed low frequency, variable frequency, 
and fixed high frequency components oriented along the X axis and referred 
to in Figure 10 by the function h(y) . The mask 3 also has another pattern 
component oriented along the line Y = Z and referred to by the function 
g(y+z) which line is designated by the reference numeral 16 and is seen to 
extend diagonally across the mask 3. That component has fixed low and 
high frequency components. Additionally, a variable frequency component 
extends along the Z axis and is designated by the reference numeral 18 as 
well as by the function j(y,z). That component comprises a variable 
frequency component . 

For this embodiment, the mask transmissivity function is: 

Mask(y,z) = g(y+z) + h(y) + j(y,z) (22) 
Replacing g(y+z), h(y), and j(y/Z) by their constituent functions, the 
mask transmissivity is described by: 

cosfr ,Z.2^| + cos[f . — •2x.\ l+- — <A/ »modl || + cos| F l .~.2/r) 
Mask( y> 7) = 0.5 + 1 — - ^ ^ - ^ =^ 

\J F .i^. 2 lJ f .^.^.Li^.^^VcJ^ t <>J2U*j| 

I I N ) \ v Wlrt W<2 N { UN )) "Xhfihi n )\ 
+ k _ (23) 

The calculations of image shift are: 

Y m _vari = [Fvariableldet*!^ - Fmidl^J 

/[Fmaxl^- Fminl^J * [3*N] pixels (24) 



A m var2 ~ 

[Fvariable2 det *K m - Fmid2 mask ] 
/[Fmax2 mask - Fmin2 roask ] * [3*N] pixels (25) 

To compute the shift of the low and high frequency components, 

corresponding to medium and fine shifts, and the function g(y+z), the 
phase at the array midpoint must be computed. The FFT routine returns the 

phase of frequencies at 0°, so to compute the midpoint phase, (i = 511.5) 



the following expressions are used: 

Arg(Flowl det ) | is5U . 5 = Arg (Flowl det ) | io0 + Mod 2n (511 . 5*Flowl det *2n) (26) 

Arg(Fhighl det ) | iB511 . 5 = Arg (Fhighl det ) | io0 + Mod 2n (511 . 5*Fhighl det *2n) (27) 

Arg(Flow2 det ) | i=511 . 5 = Arg (Flow2 det ) | i=0 + Mod 2n (511 . 5*Flow2 det *2n) (28) 

Arg(Fhigh2 det ) | i=511 . 5 = Arg (Fhigh2 det ) | ia0 + Mod 2n (511 . 5*Fhigh2 det *2n) (29) 

The shifts relative to mask functions g(y+z) and h(y-z) are 
indicated by the fixed frequency components by: 

Y m _iowi = [Plow! + (Arg(Flowl det ) | io511 . 5 /2n)] /Flowl^ (30) 

V m_highl = t^highl + 

(Arg(Fhighl det ) | i=511 . 5 /2n) ] /Fhighl^ (31) 

low2 = [Plow2 

(Arg(Flow2 det ) | io511 . 5 /2n)] /Flow2 mask (32) 

V m_high2 = t^high + 

(Arg(Fhigh2 det ) | i . 5U . 5 /2n) ] /Fhigh2 mask (33) 

Zm_hi9hl =0 (34) 

Z m _high2 = Yra_high2 (3 5) 

Where P lowl and P low2 are determined from the variable frequency shift 
estimations Y m _ varl and Y rava r2/ Phighi and P hi g h2 are determined from the low 
frequency shift estimations Y mlowl and Y mlow2 , as follows: 

P lowl = Nearest_Integer [ (Y^^^Flowl^ek) - (Arg (Flowl det ) | i=511 . 5 /2n) ] (36) 

Phi 9 h2 = Nearest_Integer [ (Y nuaow2 *Fhigh2 mask ) - (Arg (Fhigh2 det ) | i=511 . 5 /2n) ] (37) 

The overall coordinates of incidence to the mask are computed by: 

Y m = Y mhighl (3 8) 

z m = Y mhigh2 + 2Y m _ highl (3 9) 
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The angles of incidence to the mask are computed using: 



9 m = Tan" 1 (Y m /X m ) (40) 

(41) 



^yra 

8 2m = Tan' 1 (ZjXj 



As such, an invention has been disclosed in terms of preferred 
embodiments thereof which fulfill each and every one of the objects of the 
invention as set forth hereinabove, and provide a new and useful sensor 
for determining the angular position of a radiating point source in two 
dimensions, of great novelty and utility, along with the method for 
determining the angular position. 

Of course, various changes, modifications and alterations in the 
teachings of the present invention may be contemplated by those skilled in 
the art without departing from the intended spirit and scope thereof. 

As such, it is intended that the present invention only be limited 
by the terms of the appended claims. 
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